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Abstract 

The CORSIKA program, usually used to simulate extensive cosmic ray air showers, has 
been adapted to work in a water or ice medium. The adapted CORSIKA code was used 
to simulate hadronic showers produced by neutrino interactions. The simulated showers 
have been used to study the spatial distribution of the deposited energy in the showers. 
This allows a more precise determination of the acoustic signals produced by ultra high 
energy neutrinos than has been possible previously. The properties of the acoustic signals 
generated by such showers are described. 
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1 Introduction 



In recent years interest has grown in the detection of very high energy cosmic ray neutrinos [1]. 
Such particles could be produced in the cosmic particle accelerators which make the charged 
primaries or they could be produced by the interactions of the primaries with the Cosmic Mi- 
crowave Background, the so called GZK effect [2]. The flux of neutrinos expected from these 
two sources has been calculated [3,4]. It is found to be very low so that large targets are needed 
for a measurable detection rate. It is interesting to measure this neutrino flux to see if it is 
compatible with the values expected from these sources, incompatibility implying new physics. 

Searches for cosmic ray neutrinos are ongoing in AMANDA [5], IceCube [6], ANTARES 
[7] and NESTOR [8], detecting upward going muons from the Cherenkov light in either ice or 
water. In general, these experiments are sensitive to lower energies than discussed here since the 
Earth becomes opaque to neutrinos at very high energies. The experiments could detect almost 
horizontal higher energy neutrinos but have limited target volume due to the attenuation of the 
light signal in the ice. The Pierre Auger Observatory, an extended air shower array detector, will 
also search for upward and almost horizontal showers from neutrino interactions [9]. In addition 
to these detectors there are ongoing experiments to detect the neutrino interactions by either 
radio or acoustic emissions from the resulting particle showers [1]. These latter techniques, 
with much longer attenuation lengths, allow very large target volumes utilising either large 
ice fields or dry salt domes for radio or ice fields, salt domes and the oceans for the acoustic 
technique. 

In order to assess the feasibility of each technique the production of the particle shower from 
neutrino interactions needs to be simulated. Since experimental data on the interactions of such 
high energy particles do not exist it is necessary to use theoretical models to simulate them. 
The most extensive ultra high energy simulation program which has so far been developed is 
CORSIKA [10]. However, this program has been used previously only for the simulation of 
cosmic ray air showers. The program is readily available [10]. 

Different simulations are necessary for the radio and acoustic techniques. Radio emission 
occurs due to coherent Cherenkov radiation from the particles in the shower, the Askaryan 
Effect [11]. The emitted energy is sensitive to the distribution of the electron-positron asymme- 
try which develops in the shower and which grows for lower energy electromagnetic particles. 
Hence, to simulate radio emission, the electromagnetic component of the shower must be fol- 
lowed down to very low kinetic energies (~ 100 keV) [12]. In contrast, an acoustic signal 
is generated by the sudden local heating of the surrounding medium induced by the particle 
shower [13]. Thus to simulate the acoustic signal the spatial distribution of the deposited en- 
ergy is needed. Once the electromagnetic energy in the shower reaches the MeV level (electron 
range ~ 1 cm) the energy can be simply added to the total deposited energy and the simula- 
tion of such particles discontinued. Extensive simulations have been carried out for the radio 
technique [14]. However, the simulations for the acoustic technique are less advanced. Some 
work has been done [15, 16] using the Geant4 package [17]. However, this work is restricted to 
energies less than 10 5 GeV for hadron showers since the range of validity of the physics models 
in this package does not extend to higher energy hadrons. 

In this paper the energy distributions of showers produced by neutrino interactions in sea 
water at energies up to 10 12 GeV are discussed. The distributions are generated using the air 
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shower program CORSIKA [10] modified to work in a sea water medium. The salt compo- 
nent of the sea water has a negligible effect and the results are presented in distance units of 
g cm -2 , hence they should be applicable to ice also. The computed distributions have been pa- 
rameterised and this parameterisation is used to develop a simple program to simulate neutrino 
interactions and the resulting particle showers. The properties of the acoustic signals from the 
generated showers are also presented. 



2 Adaptation of the CORSIKA program to a water medium 

The air shower program, CORSIKA (version 6204) [10], has been adapted to run in sea water 
i.e. a medium of constant density of 1.025 g cm~ 3 rather than the variable density needed for 
an air atmosphere. Sea water was assumed to consist of a medium in which 66.2% of the atoms 
are hydrogen, 33.1% of the atoms are oxygen and 0.7% of the atoms are made of common salt, 
NaCl. The salt was assumed to be a material with atomic weight and atomic number A=29.2 and 
Z=14, the mean of sodium and chlorine. The purpose of this is to maintain the structure of the 
program as closely as possible to the air shower version which had two principal atmospheric 
components (oxygen and nitrogen) with a trace of argon. The presence of the salt component 
had an almost undetectable effect on the behaviour of the showers. 

Other changes made to the program to accommodate the water medium include a modifica- 
tion of the stopping power formula to allow for the density effect in water H. This only affects 
the energy loss for hadrons since the stopping powers for electrons are part of the EGS [18] 
package which is used by CORSIKA to simulate the propagation of the electromagnetic com- 
ponent of the shower. Smaller radial binning of the shower was also required since shower radii 
in water are much smaller than those in air. In addition the initial state energy for electrons 
and photons above which the LPM effect [19] was simulated in the program was reduced to the 
much lower value necessary for water[E The LPM effect suppresses pair production from high 
energy photons and bremsstrahlung from high energy electrons. Similarly, the interactions of 
neutral pions had to be simulated at lower energy than in air because of the higher density water 
medium. In all about 100 detailed changes needed to be made to the CORSIKA program to 
accommodate the water medium. 

To test the implementation of the LPM effect [19] in the program 100 showers from incident 
gamma ray photons at several different energies were generated and the mean depth of the first 
interaction (the mean free path) calculated. The observed mean free path was found to be 
in agreement with the expected behaviour when both the suppression of pair production and 
photonuclear interactions were taken into account (see Figured]). This showed that the LPM 
effect had been properly implemented in CORSIKA. 

Considerable fluctuations between showers occurred. These are expressed in terms of the 
ratio of the root mean square (RMS) deviation of a given parameter to its mean value: the 

'The shower maximum was observed to peak at a depth 2.4 ± 1.1% less in sea water than in fresh water with 
the same peak energy deposited, for protons of energy 10 5 GeV. 

2 The stopping power was computed using the Bethe-Bloch formula [20] and the density effect from the formu- 
lae of Sternheimer et al [21]. 

3 The level was set at 1 TeV compared to the characteristic energy for water Elpm = 270 TeV [20]. 
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RMS peak energy deposit to the mean peak energy deposit was observed to be 14% at 10 5 GeV 
reducing to 4% at 10 11 GeV, that for the depth of the peak position varied from 19% to 7.4% 
and for the full width at half maximum of the shower from 63% to 18%. To smooth out such 
fluctuations averages of 100 generated showers will be taken in the following. The statistical 
error on the averages is then given by these RMS values divided by 10. The hadronic energy 
contributes only about 10% to the shower energy at the shower peak, the remainder being carried 
by the electromagnetic part of the shower. 

The simulations were all carried out in a vertical column of sea water 20 m long. The 
deposited energy generated by CORSIKA was binned into 20 g cm -2 slices longitudinally 
and 1.025 g cm~ 2 annular cylinders radially for < r < 10.25 g cm" 2 and 10.25 g cm -2 
for 10.25 < r < 112.75 g cm~ 2 where r is the distance from the vertical axis. To reduce 
computing times, the thinning option was used i.e. below a certain fraction of the primary 
energy (in this case 10~ 4 ) only one of the particles emerging from the interaction is followed 
and an appropriated weight is given to it [22] . The simulation of particles continued down to cut- 
off energies of 3 MeV for electromagnetic particles and 0.3 GeV for hadrons. When a particle 
reached this cut-off, the energy was added to the slice where this occurred. The QGSJET [23] 
model was used to simulate the hadronic interactions. 



3 Comparison with other simulations 

3.1 Comparison with Geant4 

Proton showers were generated in sea water using the program Geant4 (version 8.0) [17] and 
compared with those generated in CORSIKA. Unfortunately, the range of validity of Geant4 
physics models for hadronic interactions does not extend beyond an energy of 10 5 GeV. Hence 
the comparison is restricted to energies below this. 

Figure [2] shows the longitudinal distributions of proton showers at energies of 10 4 and 10 5 
GeV (averaged over 100 showers) as determined from Geant4 and CORSIKA. The showers 
from CORSIKA tend to be slightly broader and with a smaller peak energy than those generated 
by Geant4. The difference in the peak height is ~ 5% at 10 4 GeV rising to ~ 10% at energy 10 5 
GeV. Figure [3] shows the radial distributions. The differences in the longitudinal distributions 
are reflected in the radial distributions. However, the shapes of the radial distributions are very 
similar between Geant4 and CORSIKA, with CORSIKA producing ~ 10% more energy near 
the shower axis at depths between 450 and 850 g cm~ 2 where most of the energy is deposited. 
The acoustic signal from a shower is most sensitive to the radial distribution, particularly near 
the axis (r ~ 0). It is relatively insensitive to the shape of the longitudinal distribution. 

3.2 Comparison with the simulation of Alvarez-Muniz and Zas 

The CORSIKA simulation was also compared with the longitudinal shower profile for protons 
computed in the simulation by Alvarez-Muniz and Zas (AZ) [24]. There was a reasonable 
agreement between the longitudinal shower shapes from CORSIKA and those shown in Figure 
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2 of ref. [24]. However, the numbers of electrons and positrons at the peak of the CORSIKA 
showers was ~ 20% lower than those from ref. [24]. This number is sensitive to the energy 
below which these particles are counted and this is not specified in [24]. Hence the agreement 
between CORSIKA and their simulation is probably satisfactory within this uncertainty. 

In conclusion, the modifications made to CORSIKA to simulate high energy showers in a 
water medium give results which are compatible with the predictions from the Geant4 simula- 
tions for energy less than 10 5 GeV and the simulation of AZ within 20%. This is taken to be 
the accuracy of the simulation program assuming that there are no unexpected and unknown 
interactions between the centre of mass energy explored at current accelerators and those stud- 
ied in these simulations. Studies of the sensitivity of the CORSIKA simulation to the different 
models of the hadronic interactions have been reported in reference [25]. They find that the 
peak number of electrons plus positrons varies by ~ 20% for proton showers in air depending 
on the choice of the hadron interaction model used. These differences are similar in magnitude 
to the differences between the AZ, Geant4 and CORSIKA simulations reported here. Hence 
the observed differences between the Geant4, AZ and CORSIKA simulations in water could be 
within the uncertainties of the hadronic interaction models. 

4 Simulation of neutrino induced showers 

Neutrinos interact with the nuclei of the detection medium by either the exchange of a charged 
vector boson (W + ), i.e. charged current (CC) interactions or the exchange of the neutral vector 
boson (Z°), i.e. neutral current (NC) deep inelastic scattering interactions (see for example 
[26]). The ratio of the CC to NC interaction cross sections is approximately 2:1. The CC 
interactions produce charged secondary scattered leptons while the NC interactions produce 
neutrinos. The hadron shower carries a fraction y of the energy of the incident neutrino and 
the scattered lepton the remaining fraction 1 — y. We assume that the neutrino flavours are 
homogeneously mixed when they arrive at the Earth by neutrino oscillations. Hence in the 
CC interactions electrons, // and r leptons will be produced as the scattered leptons in equal 
proportions. At the energies we shall consider, these particles behave in a manner similar to 
minimum ionising particles for fi and r leptons. This is almost true also for electrons for 
which the bremsstrahlung process will be suppressed by the LPM effect. Hence the charged 
scattered leptons contribute little to the energy producing an acoustic signal. In the case of NC 
interactions there is no contribution to this energy from the scattered lepton. For these reasons 
the contribution of the scattered lepton to the shower profile is ignored beyond z = 20 m in 
what follows. 

It is interesting to note that a r lepton can decay to hadrons or a very high energy electron or 
muon can produce bremsstrahlung photons at large distances from the interaction point. These 
can initiate further distant showers, the so called "double bang" effect. The stochastic nature of 
such electron showers is studied in [15, 16]. These effects are not considered in this study. 

4.1 Neutrino-nucleon interaction cross sections. 

A number of groups have computed the high energy neutrino-nucleon interaction cross sections, 
a, [27-29]. In the quark parton model of the nucleon for the single vector boson exchange pro- 
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cess, the differential cross section for CC interactions can be expressed in terms of the measured 
structure functions of the target nucleon F 2 and xF 3 as 

where is the Fermi weak coupling, M\y is the mass of the weak vector boson, Q 2 is the 
square of the four momentum transferred to the target nucleon, y = v jE where v is the energy 
transferred to the nucleon (y = E — E' with E and E' the energies of the incident and scattered 
leptons) and x = Q 2 /2Mv is the fraction of the momentum of the target nucleon carried by the 
struck quark (here x and y are defined for a stationary target nucleon). The plus (minus) sign 
is for neutrino (anti-neutrino) interactions. It can be seen that y is the fraction of the neutrino's 
energy which is converted into the energy of the hadron shower. A similar expression can 
be written down for the NC interaction (see for example [26]) which has a ratio to the CC 
cross section varying from 0.33 to 0.41 as the neutrino energy increases from 10 4 to 10 13 GeV. 
The structure functions F 2 and xF 3 are the sum of the quark distribution functions which have 
been parameterised by fitting data [30,31]. It can be shown that Q 2 = sxy where s = 2ME 
is the square of the centre of mass energy (M is the target nucleon mass). To compute the 
cross sections the structure functions must be calculated at values of x < Myy / s i.e. at values 
well outside the region of the fits to the parton distribution functions (PDFs) which have been 
performed for x > 10~ 5 , the range of current measurements. The extrapolation outside the 
measurement range is discussed in [27], [29] and [32,33]. Here we adopt the procedure of 
extrapolating linearly on a log-log scale from the parameterised parton distribution functions 
of [30] computed at x = 10~ 4 and x = 10 -5 . By considering various theoretical evolution 
procedures it is estimated in [29] that the procedure has an accuracy of ~ 32% per decade 
and we use this as an estimate of the accuracy of the calculation. However, this could be an 
underestimate [34]. 

The expression in equation \T\ for charged current interactions and the one for neutral current 
interactions were integrated to obtain the total neutrino-nucleon interaction cross section, the 
value of the fraction of events per interval of y, 1/ada/dy, and the mean value of y. The total 
cross section was found to be in good agreement with the values in [27, 29] and in reasonable 
agreement with [28] which is based on a model different from the quark parton model. Fig- 
ure H] shows the mean value of y obtained from this procedure (solid curve) and the effect of 
multiplying or dividing the PDFs by a factor 1.32 per decade (dashed curves) as an indication 
of the possible range of uncertainties in the extrapolation of the PDFs. Figure \5\ shows the y 
dependence of the cross section for different neutrino energies. 

4.2 A simple generator for neutrino interactions. 

A simple generator for neutrino interactions in a column of water of thickness 20 m was con- 
structed as follows. The neutrino interacts at the top of the water column (z=0, with the z axis 
along the axis of the column). The energy fraction transferred, y, for the interaction was gener- 
ated, distributed according to the curve for the energy of the neutrino shown in Figure [51 This 
allows the energy of the hadron shower to be calculated for the event. The assumption was 
made that these hadron showers will have approximately the same distributions as those of a 
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proton interaction at z=0 (see Section [4731 for a test of this assumption). A series of files of 
100 such proton interactions were generated at energies in steps of half an order of magnitude 
between 10 5 and 10 12 GeV. The hadron shower for each neutrino interaction was selected at 
random from the 100 showers in the file at the proton energy closest to the energy of the hadron 
shower. The deposited energy in each bin was then multiplied by the ratio of the energy of the 
hadron shower to that of the proton shower. This is made possible because the shower shapes 
vary slowly with shower energy. For example, the ratio of the peak energy deposit per 20 g 
cm -2 slice to the shower energy varies from 0.037 to 0.030 as the proton shower energy varies 
from 10 5 to 10 12 GeV. 



4.3 The HERWIG neutrino generator. 

The CORSIKA program has an option to simulate the interactions of neutrinos at a fixed 
point [35]. The first interaction is generated by the HERWIG package [36]. This option was 
adapted to our version of CORSIKA in sea water. Some problems were encountered with the y 
dependence of the resulting interactions due to the extrapolation of the PDFs to very small x at 
high energies. This only affects the rate of the production of the showers at different y and the 
distribution of the hadrons produced in the interaction at a given y should be unaffected. 

A total of 700 neutrino interactions were generated at an incident neutrino energy of 2 • 10 11 
GeV. These were divided into the shower energy intervals 0.5 — 2 ■ 10 10 , 2 — 4 ■ 10 10 , 4 — 7.5 • 10 10 , 
0.75 — 1.3 ■ 10 11 and 1.3 — 2 • 10 11 . The showers in which the scattered lepton energy disagreed 
with the shower energy by more than 20% were eliminated leading to a loss of 17% of the 
events with shower energy greater than 0.5 • 10 10 GeV. This is due to radiative effects and 
misidentification of the scattered lepton. Approximately 70 events remained in each energy 
interval. The energy depositions from these were averaged and compared to the averages from 
proton showers scaled by the ratio of the shower energy to the proton energy. Figure [6] shows 
the longitudinal distributions of the hadronic shower energy deposited for the different energy 
intervals (labelled E w ) compared to the scaled proton distributions. Figure [7] shows a sample 
of the transverse distributions. 

There is a good consistency between the proton and neutrino induced showers. The proton 
showers peak, on average, 20 g cm~ 2 shallower in depth with a peak energy 2% larger than the 
neutrino induced showers. This is small compared to the overall uncertainty. The slight shift in 
the longitudinal distribution is reflected as a normalisation shift in the radial distributions. We 
conclude therefore that to equate a proton induced shower starting at the neutrino interaction 
point to that from a neutrino is a satisfactory approximation. 



5 Parameterisation of showers 

In this section a parameterisation of the energy deposited by the showers generated by COR- 
SIKA (averaged over 100 showers depositing the same total energy) is described. Other avail- 
able parameterisations will then be compared with the showers generated by CORSIKA. 
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The acoustic signal generated by a hadron shower depends mainly on the energy deposited 
in the inner core of the shower. This is illustrated in figure [8] which shows the contribution to 
the acoustic signal from cores of different radii. This figure shows that it is crucial to represent 
the deposited energy well at radius less than 2.05 g cm~ 2 . The calculation of the acoustic signal 
from the deposited energy is described in section [6l 



5.1 Parameterisation of the CORSIKA Showers 

The differential energy deposited was parameterised as follows 

d 2 E 



drdz 



L(z, E L ) ■ R(r, z, E L ) (2) 



where the function L(z, El) represents the longitudinal distribution of deposited energy and 
R(r, z, El) the radial distribution. Here E L is log 10 E with E the total shower energy. 

The function L(z, E L ) = dE/dz is a modified^ version of the Gaisser-Hillas function [37]. 
This function represents the longitudinal distribution of the energy deposited. 

L(z, E L ) = P 1L ( Z ~ P2L ) "^' + ^ exp ( -) (3) 

Here the parameters P nL were fitted to quadratic functions of E L = log 10 E with values 
p 

— = 2.760 • 10~ 3 - 1.974 • 1(T 4 £ L + 7.450 • 10~ 6 E 2 L (4) 
E 

P 2L = -210.9 - 6.968 ■ 10~ 3 £ L + 0.1551/^ (5) 

P 3L = -41.50 + 113.9£ L - 4.103^2 (6) 

P 4 l = 8.012 + ll.UE L - 0.5434£f (7) 

P 5L = 0.7999 ■ 10~ 5 - 0.004843£ L + 0.0002552/^ (8) 

P 6L = 4.563 • 10~ 5 - 3.504 • 10~ 6 £ L + 1, 315 • W~ 7 E 2 L . (9) 

The parameter P 1L represents the peak energy deposited and P iL the depth in the z coordinate 
at this peak while P 2 l, Pil, Psl and P 6 l are related to the shower width and shape in z. 

The radial distribution was represented by the NKG function [37] 

R(r,z,E L ) = if (^-"(l + f-)"*^) (10) 

1 V "lR "lR J 

where the integral 



o \ Pir P\r J r(4.5 — P 2 R 



4 The modification is to replace the shape parameter A in equation 3.5 of reference [37] by the quadratic 
expression in z in equation[3] 
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The parameter P 1R was found to vary strongly with depth while P 2R was only a weak function 
of depth. The parameters P nR (with n = 1,2) were each represented by the quadratic form 



P nR = A + Bz + Cz 2 (11) 

and the quantities A, B, C parameterised as quadratic functions of E L , This gave for P 1R 

A = 0.01287£| - 0.2573£ L + 0.9636 (12) 

B = -0.4697 • 10" 4 ££ + 0.0008072£ L + 0.0005404 (13) 

C = 0.7344 • 10~ 7 ££ - 1.375 • 10~ 6 E L + 4.488 • 10" 6 (14) 
and for the parameter P 2 r 

A = -0.8905 • 10~ 3 E 2 L + 0.007727£l + 1.969 (15) 

B = 0.1173 • lO" 4 ^ - 0.0001782£ L - 5.093 • 10~ 6 (16) 

C = -0.1058 • 1Q- 7 E 2 L + 0.1524 • 10~ 6 E L - 0.1069 • 10~ 8 . (17) 



The fit was made in a depth range where dE/dz was greater than 10% of the peak value 
defined by equation|U The program MINUIT [38] was used to minimise the squared fractional 
deviations 

where F{ and Di refer to the fitted value and the value observed in the ith bin from the COR- 
SIKA showers, respectively. In order to improve the fit at small radii the contributions to x 2 
were arbitrarily weighted by 10 for r < 2.05 g cm~ 2 , 4 for 2.05 < r < 3.075 g cm~ 2 , unity 
for 3.075 < r < 51.25 g cm" 2 and 0.25 for r > 51.25 g cm" 2 . The RMS value of the frac- 
tional deviations was 3.4% for radii less than 51.25 g cm~ 2 and for energies greater than 10 6,5 
GeV. The fit becomes poorer at lower energies and greater radii than these. Integrating the pa- 
rameterisation shows that the fraction of the total energy computed from the fit within the fit 
range was 91% averaged over the deposited energy range 10 7 to 10 12 GeV. The corresponding 
fraction directly from the CORSIKA distributions was 92.5%, averaged over the same energy 
range. When applying this parameterisation at depths with smaller energy deposit than 10% of 
the peak value, the energy was assumed to be confined to an annular radius of 1.025 g cm~ 2 . 
There was a good agreement (within 5% at the peak) between the acoustic signal computed 
using this parameterisation and that taken directly from the CORSIKA showers. 



5.2 The parameterisation used by the SAUND Collaboration 

The SAUND Collaboration [39] uses the following parameterisation [40], based on the NKG 
formulae (e.g. see reference [37]), for the energy deposited per unit depth, z, and per unit 
annular thickness at radius r from a shower of energy E 

d 2 E z 

— — = Ek( )* exp (t - z/X) 2nrp{r) (19) 

(ITCIZ Zmax 
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where z max = 0.9X \n(E/E c ) is the maximum shower depth, X = 36.1 g cm 2 is the radia- 
tion length and E c = 0.0838 GeV. The constants t = z max / A where A = 130-5 log 10 (£/10 4 GeV) 
g cm -2 and k = t l ~ l / exp (t)\T(t). The radial density is given by 

r(r)=^ + a r* \J { ^; ] _ 2s) (20) 

where a = r/r M with r M = 9.04 g cm -2 , the Moliere radius in water, and s = 1.25. Figure 
[9] shows the radial distributions from CORSIKA compared with the absolute predictions of this 
p arameteri s ation . 

There is qualitative agreement between the parameterisation and the CORSIKA results. The 
difference in normalisation is explained by the somewhat different longitudinal profiles of the 
CORSIKA showers from the SAUND parameterisation. The latter are broader with a lower 
peak energy deposit and a depth of the maximum which is larger than the CORSIKA showers. 
CORSIKA predicts more energy at small r than the SAUND parameterisation. Quantitatively, 
51% of the shower energy is contained within a cylinder of radius 4 cm for the CORSIKA 
showers compared to 35% from the SAUND parameterisation. These fractions are approxi- 
mately independent of energy. Hence, in acoustic detectors a harder frequency spectrum for the 
acoustic signals is predicted by CORSIKA than by the SAUND parameterisation. Note that in 
the fit described in Section I5T1 the values of the parameter Pm (equivalent to Rm in equation 
[20]) were strongly depth dependent and much lower than the Moliere radius in water, assumed 
by the SAUND collaboration. In addition, the value of P 2 r (equivalent to s in equation l20l) 
while relatively constant tended to be at a higher value (~ 1.9) than that assumed by SAUND. 



5.3 The parameterisation used by Niess and Bertin 

Hadron showers, generated by Geant4 (version 4.06 p03), were studied up to energies of 10 5 
GeV and electromagnetic showers to higher energies by Niess and Bertin [15,16]. The hadronic 
showers were parameterised as follows. 

^ = rf{z)g{r,z) (21) 

with 

x a r» 

where E is the energy of the hadron shower, X is the radiation length in water, z' = z/X , 
b = 0.56 as determined from the fit and a is chosen to satisfy z' max = (a — l)/b. Here z' max is 
the depth in radiation lengths at which the shower maximum occurs. This is parameterised as 

4^ = 0.65 log(|4 + 3.93 (23) 
Ec 

with E c = 0.05427 GeV. The radial distribution function is parameterised as 

g(r,z)=g (^ S j (24) 
9 



where = 3.5 cm, n = n\ = 1.66 — 0.29(z/z max ) for r < r, and n = n 2 = 2.7 for r > r^. 
The constant g is chosen to be (2 — ni) (n 2 — 2)/ ((n 2 — ^i)?" 2 ) so that the integral of the radial 
distribution is unity. 

Figure [10] shows the radial distributions from this parameterisation compared with the pre- 
dictions of CORSIKA. There is quite good agreement between the two. There is a difference in 
the normalisation with depth since Geant4, on which this parameterisation is based, produces 
showers which tend to develop more slowly with depth than those from CORSIKA (see Fig- 
ure |2l). Furthermore, both this and the SAUND parameterisation (Section [531) assume a linear 
variation of the shower peak depth with log E whereas CORSIKA gives a clear parabolic shape 
(see equation[6]). This is illustrated in FigureQT] The Niess-Bertin parameterisation predicts that 
56% of the shower energy is contained within a cylinder of radius 4 cm in reasonable agreement 
with the value of 51% from CORSIKA (these values are almost independent of energy). 



6 The acoustic signals from the showers. 

The pressure, P, from a hadron shower depositing total energy E at time t resulting from the 
deposition of relative energy density e = (l/E)(l/2irr)d 2 E/drdz at a point distant d from the 
volume, dV, follows the form [13] 

p ^-£-JiJ t { Ht - i/c) ) iV <25) 

where the integral is over the total volume of the shower. Here (3 = 2.0 ■ 10~ 4 is the thermal 
expansion coefficient of the medium at 14°C, C p = 3.8 • 10 3 J kg -1 K 1 is the specific heat 
capacity and c = 1500 ms _1 is the velocity of sound in the sea water. 

Acoustic signals seen by an observer at distance r from the shower centre are computed from 
equation (1251) as follows. Points are produced randomly throughout the volume of the shower 
with density proportional to the deposited energy density and the time of flight from every 
produced point to the observer calculated. The flight times to the observer are histogrammed 
over 2 n bins (in this case n = 10 is chosen) centred on the mean flight time and with a suitable 
bin width, r (chosen here to be I (is). The counts in each bin of the histogram are divided by r 
yielding the function E xyz (t). The Fourier transform of the pressure wave is then 

1 r°° EQ d 1 E B r°° 1 E Q 

(26) 

using the standard Fourier transform theorem, that taking the derivative in the time domain is the 
same as multiplying by icu in the frequency domain. The Fourier transform E xyz {uS) at angular 
frequency uj is evaluated numerically by a fast Fourier Transform (FFT) from the histogram 
E xy z(t)- A correction is applied for attenuation in the water by a factor A(u) = e~ a ^ r where 
a(uj) is the frequency dependent attenuation coefficient. The pressure as a function of time is 
then evaluated numerically by an inverse FFT using frequency steps from zero to the sampling 
frequency (the inverse of the bin width r i.e. 1 MHz in this case). This gives 

1 n=511 

jP(t) = 1024 £ P M A Me inQ (27) 

n=-512 
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where Q = 27r/1024 radians and uj n /2n = nQ/2n MHz is the nth frequency. The attenuation 
coefficient a{uS) is computed either according to the formulae in [42] or using the complex 
attenuation given in [15, 16]. This method of calculation was computationally much faster than 
the evaluation of the space integral given in equation 18 of reference [13] and gave identical 
results. 

Acoustic pulses, computed with the complex attenuation described in [15, 16], using the 
parameterisations of the shower profile given above are shown in Figure [[21 It can be seen that 
the parameterisation developed here gives similar results to that described in [15, 16] despite 
the fact that the latter was an extrapolation from low energy simulations. The parameterisation 
used by SAUND [39,40] gives smaller signals concentrated at somewhat lower frequencies. 

Further properties of the acoustic signals are shown in Figures [13] to [161 The pulses tend 
to be somewhat asymmetric with the asymmetry defined by |P maa; | — \P m in\ / \P ma x\ + \Pmin\- 
The complex nature of the attenuation enhances this asymmetry. This is most evident in the far 
field conditions e.g. at 5km where non complex attenuation would yield a totally symmetric 
pulse. Figure [13] shows the angular dependence of the peak pressure. Here the angle is that 
subtended by the acoustic detector relative to the plane, termed the median plane, through the 
shower maximum at right angles to the axis of the shower. The parameterisation derived here 
gives a somewhat narrower angular spread than the others. This could be due to the slightly 
longer showers predicted by CORSIKA than the others. Figure [13] also shows the asymmetry of 
the pulse as a function of this angle. The pulse initially becomes more symmetric moving out of 
the median plane and then the asymmetry becomes negative at larger angles. Figure [T41 shows 
the decrease of the pulsed peak pressure with distance from the shower in the median plane and 
the asymmetry with distance in this plane. Figures [15] and [16] show the frequency composition 
of the pulses at different angles to the median plane at 1 km from the shower and at different 
distance in the median plane, respectively. 

7 Conclusions 

The simulation program for high energy cosmic ray air showers, CORSIKA, has been modified 
to work in a water or ice medium. This allows both hadron and neutrino showers to be generated 
in the medium over a wide range of energy (10 5 to 10 12 GeV). The properties of hadronic 
showers in water simulated by CORSIKA agree with those from other simulations to within 
10 — 20%. A similar uncertainty has been noted previously from the variations in CORSIKA 
showers in air generated by different models of the hadron interactions. However, none of 
the other available simulations for water cover the range of energies accessible to CORSIKA. 
The hadronic showers produced by neutrino interactions are shown to have similar profiles to 
proton showers which deposit the same amount of energy to that from the neutrino and which 
start at the interaction point of the neutrino. The properties of the neutrino interactions are 
described. A parameterisation of the shower profiles generated by CORSIKA is given. There is 
reasonable agreement with the parameterisation based on the Geant4 simulations at low energy 
(< 10 5 GeV) developed by Niess and Bertin. However, the agreement with the parameterisation 
used by the SAUND Collaboration, which is based on the NKG formalism, is less good. The 
position of the shower maximum, determined from the CORSIKA program, is found to vary 
quadratically with \ogE rather than linearly as assumed in the latter two parameterisations. 
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The acoustic signals generated by neutrino interactions using CORSIKA and by the two 
other parameterisations are described and their properties are studied. The acoustic signal is 
found to be very sensitive to the energy deposited close to the shower axis. 
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Interaction length of photons in water (CORSIKA) 
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Figure 1 : The interaction length for high energy gamma rays versus the photon energy measured 
in CORSIKA (data points with statistical errors). The dash dotted curve shows the pair produc- 
tion length computed from the LPM effect using the formulae of Migdal [19]. The solid curve 
shows the computed total interaction length, including both pair production and photonuclear 
interactions with the cross section from CORSIKA. The dashed line labelled 9/7X shows the 
expected pair production length without the LPM effect. Here X is the radiation length of the 
material. 
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Average dE/dz of 100 proton showers 
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Figure 2: Averaged longitudinal energy deposited per unit path length of 100 proton showers 
at energy 10 4 GeV (upper plot) and 10 5 GeV (lower plot) generated in Geant4 and CORSIKA 
versus depth in the water. 
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Figure 3: Averaged radial energy deposited per 20 g cm~ 2 vertical slice per unit radial distance 
for 100 proton showers at energy 10 4 GeV (left hand plots) and 10 5 GeV (right hand plots) 
generated in Geant4 and CORSIKA versus distance from the axis in the water for different 
depths of the shower. 
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Average y for different structure function extrapolations 
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Figure 4: The mean value of y as a function of energy for interactions computed according 
to the standard model with the PDFs of MRS99 [30], extrapolating x and Q 2 out of the fit 
range from x = 1CT 4 linearly on a log-log scale. The upper dashed curve shows the result of 
multiplying the PDFs by 1.32 log ( 10 / x ) for PDFs with x < 1CT 4 and the lower dashed curve by 
dividing by this factor. The deviations of the dashed curve from the solid one is an indication 
of the precision of the standard model. 
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Figure 5: The fraction of events per unit y interval for different energies computed by inte- 
grating the expressions for the CC and NC cross sections. 
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Proton (dashed) v (solid) showers of same energy 
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Figure 6: The longitudinal distribution of the deposited energy for neutrino showers (solid) 
generated by the Herwig-CORSIKA package and proton showers (dashed) scaled to the same 
values of shower energy E w . The scaling factors applied to the average of the protons showers 
with energy 10 10 GeV were 1.0 and 3.0 for E w = 10 10 GeV and E w = 3 ■ 10 10 GeV, respec- 
tively. Those applied to proton showers with energy 10 11 GeV were 0.575, 1.0 and 1.65 for 
E w = 5.75 • 10 10 E w = 10 11 and E w = 1.65 • 10 11 GeV, respectively. 
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Figure 7: The solid curves show the averaged radial energy deposited per 20 g cm~ 2 vertical 
slice per unit radial distance for 70 neutrino showers with energy transfer E w = 3 • 10 10 GeV 
(left hand plots) and E w = 1.65 • 10 11 GeV (right hand plots). The incident neutrino energy 
was 2 • 10 11 GeV. For comparison the dashed curves show the distributions from proton showers 
scaled to these energies. In the left (right) hand plots protons of energy 10 10 GeV (10 11 GeV) 
were scaled by a factor of 3 (1.65). 
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Figure 8: The acoustic signal at a distance of 1 km from the shower axis in the median plane 
computed from the average of 100 CORSIKA showers each depositing a total energy of 10 9 
GeV in the water. The dotted, dashed and solid curves shows the signals computed from the 
deposited energies within cores of radius 1.025 g cm~ 2 , 2.05 g cm -2 and the whole shower 
(solid curve), respectively. It can be seen that most of the amplitude of the signal comes from 
the energy within a core of radius 2.05 g cm~ 2 . 
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Figure 9: The radial distributions of the deposited energy at different depths from CORSIKA 
compared to the parameterisation used by the SAUND collaboration [39] for 10 6 GeV and 10 11 
GeV proton induced showers. 
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Figure 10: The radial distributions of the deposited energy at different depths from CORSIKA 
compared to the parameterisation used by the Niess and Bertin [15, 16] (labelled NB parame- 
terisation) for 10 6 GeV and 10 11 GeV proton induced showers. 
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Figure 1 1 : The depth of the shower peak as a function of log 10 E from CORSIKA (black points) 
for the showers starting at z — 0. The solid curve shows the parameterisation according to 
equation |6j The dashed (dash dotted) lines show the values assumed by the SAUND (Niess- 
Bertin) Collaborations. 
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The Acoustic Pulse from a 10 11 GeV Shower at 1km Relative Power Spectrum of a Pulse from a 10 11 GeV Shower at 1km 
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Figure 12: The left hand plot shows acoustic pulses generated from the parameterisation de- 



scribed in section I5TT1 labelled ACoRNE, the parameterisation from reference [15, 16] labelled 
NB and that from reference [39,40] labelled SAUND. These pulses were evaluated for a hadron 
shower from a neutrino interaction depositing hadronic energy of 10 11 GeV 1 km distant from 
an acoustic detector in a plane perpendicular to the shower axis at the shower maximum (the 
median plane). The right hand plot shows the frequency decomposition of the pulses in the left 
hand plot. 
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Peak Pressure vs Angle from a l0 11 GeV Shower at 1km The Asymmetry of a Pulse from a 10 GeV Shower 
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Figure 13: The left hand plot shows the variation of the peak pressure in the pulse with angle 
from the median plane at 1 km from the shower. The right hand plot shows the variation of 
the pulse asymmetry with this angle at the same distance. The curves were computed from the 
parameterisations labelled. 
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Figure 14: The left hand plot shows the decrease of the pulse peak pressure and the right hand 
plot the pulse asymmetry, both in the median plane, as a function of distance from the shower 
computed from the parameterisations. 
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Relative Power Spectrum of a Pulse from a 10 11 GeV Shower at 1km Relative Power Spectrum of a Pulse from a 10 11 GeV Shower at 1 km 
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Figure 15: The left hand plot shows the frequency decomposition of the acoustic signal, com- 
puted from the parameterisation of the CORSIKA showers, at different angles to the median 
plane at a distance of 1km from the shower and the right hand plot shows the cumulative fre- 
quency spectrum i.e. the integral of the left hand plot. 



Relative Power Spectrum of a Pulse from a 10 11 GeV Shower Relative Power Spectrum of a Pulse from a 10 11 GeV Shower 
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Figure 16: The left hand plot shows the frequency decomposition of the acoustic signal, com- 
puted from the parameterisation of the CORSIKA showers, at different distances from the hy- 
drophone in the median plane and the right hand plot shows the cumulative frequency spectrum 
i.e. the integral of the left hand plot. 
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Abstract 

The CORSIKA program, usually used to simulate extensive cosmic ray air showers, has 
been adapted to work in a water or ice medium. The adapted CORSIKA code was used 
to simulate hadronic showers produced by neutrino interactions. The simulated showers 
have been used to study the spatial distribution of the deposited energy in the showers. 
This allows a more precise determination of the acoustic signals produced by ultra high 
energy neutrinos than has been possible previously. The properties of the acoustic signals 
generated by such showers are described. 



(Submitted to Astroparticle Physics) 



1 Introduction 



In recent years interest has grown in the detection of very high energy cosmic ray neutrinos [1]. 
Such particles could be produced in the cosmic particle accelerators which make the charged 
primaries or they could be produced by the interactions of the primaries with the Cosmic Mi- 
crowave Background, the so called GZK effect [2]. The flux of neutrinos expected from these 
two sources has been calculated [3,4]. It is found to be very low so that large targets are needed 
for a measurable detection rate. It is interesting to measure this neutrino flux to see if it is 
compatible with the values expected from these sources, incompatibility implying new physics. 

Searches for cosmic ray neutrinos are ongoing in AMANDA [5], IceCube [6], ANTARES 
[7] and NESTOR [8], detecting upward going muons from the Cherenkov light in either ice or 
water. In general, these experiments are sensitive to lower energies than discussed here since the 
Earth becomes opaque to neutrinos at very high energies. The experiments could detect almost 
horizontal higher energy neutrinos but have limited target volume due to the attenuation of the 
light signal in the ice. The Pierre Auger Observatory, an extended air shower array detector, will 
also search for upward and almost horizontal showers from neutrino interactions [9]. In addition 
to these detectors there are ongoing experiments to detect the neutrino interactions by either 
radio or acoustic emissions from the resulting particle showers [1]. These latter techniques, 
with much longer attenuation lengths, allow very large target volumes utilising either large 
ice fields or dry salt domes for radio or ice fields, salt domes and the oceans for the acoustic 
technique. 

In order to assess the feasibility of each technique the production of the particle shower from 
neutrino interactions needs to be simulated. Since experimental data on the interactions of such 
high energy particles do not exist it is necessary to use theoretical models to simulate them. 
The most extensive ultra high energy simulation program which has so far been developed is 
CORSIKA [10]. However, this program has been used previously only for the simulation of 
cosmic ray air showers. The program is readily available [10]. 

Different simulations are necessary for the radio and acoustic techniques. Radio emission 
occurs due to coherent Cherenkov radiation from the particles in the shower, the Askaryan 
Effect [11]. The emitted energy is sensitive to the distribution of the electron-positron asymme- 
try which develops in the shower and which grows for lower energy electromagnetic particles. 
Hence, to simulate radio emission, the electromagnetic component of the shower must be fol- 
lowed down to very low kinetic energies (~ 100 keV) [12]. In contrast, an acoustic signal 
is generated by the sudden local heating of the surrounding medium induced by the particle 
shower [13]. Thus to simulate the acoustic signal the spatial distribution of the deposited en- 
ergy is needed. Once the electromagnetic energy in the shower reaches the MeV level (electron 
range ~ 1 cm) the energy can be simply added to the total deposited energy and the simula- 
tion of such particles discontinued. Extensive simulations have been carried out for the radio 
technique [14]. However, the simulations for the acoustic technique are less advanced. Some 
work has been done [15, 16] using the Geant4 package [17]. However, this work is restricted to 
energies less than 10 5 GeV for hadron showers since the range of validity of the physics models 
in this package does not extend to higher energy hadrons. 

In this paper the energy distributions of showers produced by neutrino interactions in sea 
water at energies up to 10 12 GeV are discussed. The distributions are generated using the air 
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shower program CORSIKA [10] modified to work in a sea water medium. The salt compo- 
nent of the sea water has a negligible effect 1 and the results are presented in distance units of 
g cm -2 , hence they should be applicable to ice also. The computed distributions have been pa- 
rameterised and this parameterisation is used to develop a simple program to simulate neutrino 
interactions and the resulting particle showers. The properties of the acoustic signals from the 
generated showers are also presented. 



2 Adaptation of the CORSIKA program to a water medium 

The air shower program, CORSIKA (version 6204) [10], has been adapted to run in sea water 
i.e. a medium of constant density of 1.025 g cm~ 3 rather than the variable density needed for 
an air atmosphere. Sea water was assumed to consist of a medium in which 66.2% of the atoms 
are hydrogen, 33.1% of the atoms are oxygen and 0.7% of the atoms are made of common salt, 
NaCl. The salt was assumed to be a material with atomic weight and atomic number A=29.2 and 
Z=14, the mean of sodium and chlorine. The purpose of this is to maintain the structure of the 
program as closely as possible to the air shower version which had two principal atmospheric 
components (oxygen and nitrogen) with a trace of argon. The presence of the salt component 
had an almost undetectable effect on the behaviour of the showers. 

Other changes made to the program to accommodate the water medium include a modifica- 
tion of the stopping power formula to allow for the density effect in water 2 . This only affects 
the energy loss for hadrons since the stopping powers for electrons are part of the EGS [18] 
package which is used by CORSIKA to simulate the propagation of the electromagnetic com- 
ponent of the shower. Smaller radial binning of the shower was also required since shower radii 
in water are much smaller than those in air. In addition the initial state energy for electrons 
and photons above which the LPM effect [19] was simulated in the program was reduced to the 
much lower value necessary for water 3 . The LPM effect suppresses pair production from high 
energy photons and bremsstrahlung from high energy electrons. Similarly, the interactions of 
neutral pions had to be simulated at lower energy than in air because of the higher density water 
medium. In all about 100 detailed changes needed to be made to the CORSIKA program to 
accommodate the water medium. 

To test the implementation of the LPM effect [19] in the program 100 showers from incident 
gamma ray photons at several different energies were generated and the mean depth of the first 
interaction (the mean free path) calculated. The observed mean free path was found to be 
in agreement with the expected behaviour when both the suppression of pair production and 
photonuclear interactions were taken into account (see Figure 1). This showed that the LPM 
effect had been properly implemented in CORSIKA. 

Considerable fluctuations between showers occurred. These are expressed in terms of the 
ratio of the root mean square (RMS) deviation of a given parameter to its mean value: the 

'The shower maximum was observed to peak at a depth 2.4 ± 1.1% less in sea water than in fresh water with 
the same peak energy deposited, for protons of energy 10 5 GeV. 

2 The stopping power was computed using the Bethe-Bloch formula [20] and the density effect from the formu- 
lae of Sternheimer et al [21]. 

3 The level was set at 1 TeV compared to the characteristic energy for water Elpm = 270 TeV [20]. 



2 



RMS peak energy deposit to the mean peak energy deposit was observed to be 14% at 10 5 GeV 
reducing to 4% at 10 11 GeV, that for the depth of the peak position varied from 19% to 7.4% 
and for the full width at half maximum of the shower from 63% to 18%. To smooth out such 
fluctuations averages of 100 generated showers will be taken in the following. The statistical 
error on the averages is then given by these RMS values divided by 10. The hadronic energy 
contributes only about 10% to the shower energy at the shower peak, the remainder being carried 
by the electromagnetic part of the shower. 

The simulations were all carried out in a vertical column of sea water 20 m long. The 
deposited energy generated by CORSIKA was binned into 20 g cm -2 slices longitudinally 
and 1.025 g cm~ 2 annular cylinders radially for < r < 10.25 g cm 2 and 10.25 g cm 2 
for 10.25 < r < 112.75 g cm~ 2 where r is the distance from the vertical axis. To reduce 
computing times, the thinning option was used i.e. below a certain fraction of the primary 
energy (in this case 10~ 4 ) only one of the particles emerging from the interaction is followed 
and an appropriated weight is given to it [22] . The simulation of particles continued down to cut- 
off energies of 3 MeV for electromagnetic particles and 0.3 GeV for hadrons. When a particle 
reached this cut-off, the energy was added to the slice where this occurred. The QGSJET [23] 
model was used to simulate the hadronic interactions. 



3 Comparison with other simulations 

3.1 Comparison with Geant4 

Proton showers were generated in sea water using the program Geant4 (version 8.0) [17] and 
compared with those generated in CORSIKA. Unfortunately, the range of validity of Geant4 
physics models for hadronic interactions does not extend beyond an energy of 10 5 GeV. Hence 
the comparison is restricted to energies below this. 

Figure 2 shows the longitudinal distributions of proton showers at energies of 10 4 and 10 5 
GeV (averaged over 100 showers) as determined from Geant4 and CORSIKA. The showers 
from CORSIKA tend to be slightly broader and with a smaller peak energy than those generated 
by Geant4. The difference in the peak height is ~ 5% at 10 4 GeV rising to ~ 10% at energy 10 5 
GeV. Figure 3 shows the radial distributions. The differences in the longitudinal distributions 
are reflected in the radial distributions. However, the shapes of the radial distributions are very 
similar between Geant4 and CORSIKA, with CORSIKA producing ~ 10% more energy near 
the shower axis at depths between 450 and 850 g cm~ 2 where most of the energy is deposited. 
The acoustic signal from a shower is most sensitive to the radial distribution, particularly near 
the axis (r ~ 0). It is relatively insensitive to the shape of the longitudinal distribution. 

3.2 Comparison with the simulation of Alvarez-Muniz and Zas 

The CORSIKA simulation was also compared with the longitudinal shower profile for protons 
computed in the simulation by Alvarez-Muniz and Zas (AZ) [24]. There was a reasonable 
agreement between the longitudinal shower shapes from CORSIKA and those shown in Figure 
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2 of ref. [24]. However, the numbers of electrons and positrons at the peak of the CORSIKA 
showers was ~ 20% lower than those from ref. [24]. This number is sensitive to the energy 
below which these particles are counted and this is not specified in [24]. Hence the agreement 
between CORSIKA and their simulation is probably satisfactory within this uncertainty. 

In conclusion, the modifications made to CORSIKA to simulate high energy showers in a 
water medium give results which are compatible with the predictions from the Geant4 simula- 
tions for energy less than 10 5 GeV and the simulation of AZ within 20%. This is taken to be 
the accuracy of the simulation program assuming that there are no unexpected and unknown 
interactions between the centre of mass energy explored at current accelerators and those stud- 
ied in these simulations. Studies of the sensitivity of the CORSIKA simulation to the different 
models of the hadronic interactions have been reported in reference [25]. They find that the 
peak number of electrons plus positrons varies by ~ 20% for proton showers in air depending 
on the choice of the hadron interaction model used. These differences are similar in magnitude 
to the differences between the AZ, Geant4 and CORSIKA simulations reported here. Hence 
the observed differences between the Geant4, AZ and CORSIKA simulations in water could be 
within the uncertainties of the hadronic interaction models. 

4 Simulation of neutrino induced showers 

Neutrinos interact with the nuclei of the detection medium by either the exchange of a charged 
vector boson (W + ), i.e. charged current (CC) interactions or the exchange of the neutral vector 
boson (Z°), i.e. neutral current (NC) deep inelastic scattering interactions (see for example 
[26]). The ratio of the CC to NC interaction cross sections is approximately 2:1. The CC 
interactions produce charged secondary scattered leptons while the NC interactions produce 
neutrinos. The hadron shower carries a fraction y of the energy of the incident neutrino and 
the scattered lepton the remaining fraction 1 — y. We assume that the neutrino flavours are 
homogeneously mixed when they arrive at the Earth by neutrino oscillations. Hence in the 
CC interactions electrons, ji and r leptons will be produced as the scattered leptons in equal 
proportions. At the energies we shall consider, these particles behave in a manner similar to 
minimum ionising particles for \i and r leptons. This is almost true also for electrons for 
which the bremsstrahlung process will be suppressed by the LPM effect. Hence the charged 
scattered leptons contribute little to the energy producing an acoustic signal. In the case of NC 
interactions there is no contribution to this energy from the scattered lepton. For these reasons 
the contribution of the scattered lepton to the shower profile is ignored beyond z = 20 m in 
what follows. 

It is interesting to note that a r lepton can decay to hadrons or a very high energy electron or 
muon can produce bremsstrahlung photons at large distances from the interaction point. These 
can initiate further distant showers, the so called "double bang" effect. The stochastic nature of 
such electron showers is studied in [15, 16]. These effects are not considered in this study. 

4.1 Neutrino-nucleon interaction cross sections. 

A number of groups have computed the high energy neutrino-nucleon interaction cross sections, 
a, [27-29]. In the quark parton model of the nucleon for the single vector boson exchange pro- 
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cess, the differential cross section for CC interactions can be expressed in terms of the measured 
structure functions of the target nucleon F 2 and xF s as 

]2 f<2 / M 2 \ 2 

WT« = ^y{l?Tk) (^,Q^-v + vV2)±v(i-vn) X F 3{ ^ )) (i) 

where Gp is the Fermi weak coupling, M w is the mass of the weak vector boson, Q 2 is the 
square of the four momentum transferred to the target nucleon, y = vjE where v is the energy 
transferred to the nucleon (v = E — E' with E and E' the energies of the incident and scattered 
leptons) and x = Q 2 /2Mv is the fraction of the momentum of the target nucleon carried by the 
struck quark (here x and y are defined for a stationary target nucleon). The plus (minus) sign 
is for neutrino (anti-neutrino) interactions. It can be seen that y is the fraction of the neutrino's 
energy which is converted into the energy of the hadron shower. A similar expression can 
be written down for the NC interaction (see for example [26]) which has a ratio to the CC 
cross section varying from 0.33 to 0.41 as the neutrino energy increases from 10 4 to 10 13 GeV. 
The structure functions F 2 and xF s are the sum of the quark distribution functions which have 
been parameterised by fitting data [30,31]. It can be shown that Q 2 = sxy where s = 2ME 
is the square of the centre of mass energy (M is the target nucleon mass). To compute the 
cross sections the structure functions must be calculated at values of x < M^/s i.e. at values 
well outside the region of the fits to the parton distribution functions (PDFs) which have been 
performed for x > 10 5 , the range of current measurements. The extrapolation outside the 
measurement range is discussed in [27], [29] and [32,33]. Here we adopt the procedure of 
extrapolating linearly on a log-log scale from the parameterised parton distribution functions 
of [30] computed at x = 10~ 4 and x = 10~ 5 . By considering various theoretical evolution 
procedures it is estimated in [29] that the procedure has an accuracy of ~ 32% per decade 
and we use this as an estimate of the accuracy of the calculation. However, this could be an 
underestimate [34]. 

The expression in equation 1 for charged current interactions and the one for neutral current 
interactions were integrated to obtain the total neutrino-nucleon interaction cross section, the 
value of the fraction of events per interval of y, 1/ada/dy, and the mean value of y. The total 
cross section was found to be in good agreement with the values in [27, 29] and in reasonable 
agreement with [28] which is based on a model different from the quark parton model. Fig- 
ure 4 shows the mean value of y obtained from this procedure (solid curve) and the effect of 
multiplying or dividing the PDFs by a factor 1.32 per decade (dashed curves) as an indication 
of the possible range of uncertainties in the extrapolation of the PDFs. Figure 5 shows the y 
dependence of the cross section for different neutrino energies. 

4.2 A simple generator for neutrino interactions. 

A simple generator for neutrino interactions in a column of water of thickness 20 m was con- 
structed as follows. The neutrino interacts at the top of the water column (z=0, with the z axis 
along the axis of the column). The energy fraction transferred, y, for the interaction was gener- 
ated, distributed according to the curve for the energy of the neutrino shown in Figure 5. This 
allows the energy of the hadron shower to be calculated for the event. The assumption was 
made that these hadron showers will have approximately the same distributions as those of a 
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proton interaction at z=0 (see Section 4.3 for a test of this assumption). A series of files of 
100 such proton interactions were generated at energies in steps of half an order of magnitude 
between 10 5 and 10 12 GeV. The hadron shower for each neutrino interaction was selected at 
random from the 100 showers in the file at the proton energy closest to the energy of the hadron 
shower. The deposited energy in each bin was then multiplied by the ratio of the energy of the 
hadron shower to that of the proton shower. This is made possible because the shower shapes 
vary slowly with shower energy. For example, the ratio of the peak energy deposit per 20 g 
cm -2 slice to the shower energy varies from 0.037 to 0.030 as the proton shower energy varies 
from 10 5 to 10 12 GeV. 

4.3 The HERWIG neutrino generator. 

The CORSIKA program has an option to simulate the interactions of neutrinos at a fixed 
point [35]. The first interaction is generated by the HERWIG package [36]. This option was 
adapted to our version of CORSIKA in sea water. Some problems were encountered with the y 
dependence of the resulting interactions due to the extrapolation of the PDFs to very small x at 
high energies. This only affects the rate of the production of the showers at different y and the 
distribution of the hadrons produced in the interaction at a given y should be unaffected. 

A total of 700 neutrino interactions were generated at an incident neutrino energy of 2 • 10 11 
GeV. These were divided into the shower energy intervals 0.5 — 2 • 10 10 , 2 — 4 • 10 10 , 4 — 7.5 • 10 10 , 
0.75 — 1.3 • 10 11 and 1.3 — 2 • 10 11 . The showers in which the scattered lepton energy disagreed 
with the shower energy by more than 20% were eliminated leading to a loss of 17% of the 
events with shower energy greater than 0.5 • 10 10 GeV. This is due to radiative effects and 
misidentification of the scattered lepton. Approximately 70 events remained in each energy 
interval. The energy depositions from these were averaged and compared to the averages from 
proton showers scaled by the ratio of the shower energy to the proton energy. Figure 6 shows 
the longitudinal distributions of the hadronic shower energy deposited for the different energy 
intervals (labelled E w ) compared to the scaled proton distributions. Figure 7 shows a sample 
of the transverse distributions. 

There is a good consistency between the proton and neutrino induced showers. The proton 
showers peak, on average, 20 g cm -2 shallower in depth with a peak energy 2% larger than the 
neutrino induced showers. This is small compared to the overall uncertainty. The slight shift in 
the longitudinal distribution is reflected as a normalisation shift in the radial distributions. We 
conclude therefore that to equate a proton induced shower starting at the neutrino interaction 
point to that from a neutrino is a satisfactory approximation. 



5 Parameterisation of showers 

In this section a parameterisation of the energy deposited by the showers generated by COR- 
SIKA (averaged over 100 showers depositing the same total energy) is described. Other avail- 
able parameterisations will then be compared with the showers generated by CORSIKA. 
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The acoustic signal generated by a hadron shower depends mainly on the energy deposited 
in the inner core of the shower. This is illustrated in figure 8 which shows the contribution to 
the acoustic signal from cores of different radii. This figure shows that it is crucial to represent 
the deposited energy well at radius less than 2.05 g cm~ 2 . The calculation of the acoustic signal 
from the deposited energy is described in section 6. 



5.1 Parameterisation of the CORSIKA Showers 

The differential energy deposited was parameterised as follows 

d 2 E 



drdz 



L(z,E L )-R(r,z,E L ) (2) 



where the function L(z, E L ) represents the longitudinal distribution of deposited energy and 
R(r, z, E L ) the radial distribution. Here E L is log 10 E with E the total shower energy. 

The function L(z, E L ) = dE/dz is a modified 4 version of the Gaisser-Hillas function [37]. 
This function represents the longitudinal distribution of the energy deposited. 

(p 3£ -p 2 l) 

Here the parameters P nL were fitted to quadratic functions of E L = log 10 E with values 
p 

= 2.760 • 10~ 3 - 1.974 • 10~ 4 £ L + 7.450 • 1(T 6 E£ (4) 

E 

P 2L = -210.9 - 6.968 • 10~ 3 E L + 0.1551#£ (5) 

P 3L = -41.50 + 113.9E L - 4.103^£ (6) 

P 4L = 8.012 + ll.UE L - 0.5434E| (7) 

p 5L = 0.7999 • 10" 5 - 0.004843£ L + 0.0002552E 2 (8) 

P 6L = 4.563 • 10" 5 - 3.504 • lO -6 ^ + 1, 315 • 10" 7 ^. (9) 

The parameter P 1L represents the peak energy deposited and P 3L the depth in the z coordinate 
at this peak while P 2L , P 4L , P 5L and P 6L are related to the shower width and shape in z. 

The radial distribution was represented by the NKG function [37] 



(10) 



where the integral 



i=i ( (i^r^ 1 + ^) (P2 - 4 - 5) ) dr = pj {4 - 5 - 2p **n p **) 



Pir' v Pir' J T(A.h-P 2R ) 



4 The modification is to replace the shape parameter A in equation 3.5 of reference [37] by the quadratic expres- 
sion in z in equation 3. 



The parameter Pi R was found to vary strongly with depth while P 2R was only a weak function 
of depth. The parameters P nR (with n = 1,2) were each represented by the quadratic form 



P nR = A + Bz + Cz 2 (11) 

and the quantities A, B>, C parameterised as quadratic functions of E L . This gave for Pi R 

A = 0.01287E| - 0.2573£ L + 0.9636 (12) 

B = -0.4697 • 10~ 4 E 2 L + 0.0008072£ L + 0.0005404 (13) 

C = 0.7344 • 10- 7 E 2 L - 1.375 • 10~ 6 E L + 4.488 • 10" 6 (14) 
and for the parameter P 2R 

A = -0.8905 • 10" 3 E| + 0.007727£ L + 1.969 (15) 

B = 0.1173 • 10~ 4 E 2 L - 0.0001782£ L - 5.093 • 10" 6 (16) 

C = -0.1058 • 1Q- 7 E 2 L + 0.1524 • 10~ 6 E L - 0.1069 • 10" 8 . (17) 



The fit was made in a depth range where dE/dz was greater than 10% of the peak value 
defined by equation 4. The program MINUIT [38] was used to minimise the squared fractional 
deviations 

where and D { refer to the fitted value and the value observed in the ith bin from the COR- 
SIKA showers, respectively. In order to improve the fit at small radii the contributions to x 2 
were arbitrarily weighted by 10 for r < 2.05 g cm" 2 , 4 for 2.05 < r < 3.075 g cm~ 2 , unity 
for 3.075 < r < 51.25 g cm" 2 and 0.25 for r > 51.25 g cm" 2 . The RMS value of the frac- 
tional deviations was 3.4% for radii less than 51.25 g cm -2 and for energies greater than 10 6 5 
GeV. The fit becomes poorer at lower energies and greater radii than these. Integrating the pa- 
rameterisation shows that the fraction of the total energy computed from the fit within the fit 
range was 91% averaged over the deposited energy range 10 7 to 10 12 GeV. The corresponding 
fraction directly from the CORSIKA distributions was 92.5%, averaged over the same energy 
range. When applying this parameterisation at depths with smaller energy deposit than 10% of 
the peak value, the energy was assumed to be confined to an annular radius of 1.025 g cm -2 . 
There was a good agreement (within 5% at the peak) between the acoustic signal computed 
using this parameterisation and that taken directly from the CORSIKA showers. 



5.2 The parameterisation used by the SAUND Collaboration 

The SAUND Collaboration [39] uses the following parameterisation [40], based on the NKG 
formulae (e.g. see reference [37]), for the energy deposited per unit depth, z, and per unit 
annular thickness at radius r from a shower of energy E 

= Ek(—Yexp(t- z/\) 2nrp(r) (19) 



where z max = 0.9X \n(E/E c ) is the maximum shower depth, X = 36.1 g cm 2 is the radia- 
tion length and E c = 0.0838 GeV. The constants t = z max /X where A = 130-5 log 10 (E/10 4 GeV) 
g cm -2 and k = t l ~ l / exp (t)XT(t). The radial density is given by 

( r ) = —a s - 2 (l + a) 8 ' 4 - 5 - fj (20) 

P[ ) r 2 / [ + ' 27rr( S )r(4.5 - 2s) ^ 

where a = r/r M with r M = 9.04 g cm~ 2 , the Moliere radius in water, and s = 1.25. Figure 
9 shows the radial distributions from CORSIKA compared with the absolute predictions of this 
parameterisation. 

There is qualitative agreement between the parameterisation and the CORSIKA results. The 
difference in normalisation is explained by the somewhat different longitudinal profiles of the 
CORSIKA showers from the SAUND parameterisation. The latter are broader with a lower 
peak energy deposit and a depth of the maximum which is larger than the CORSIKA showers. 
CORSIKA predicts more energy at small r than the SAUND parameterisation. Quantitatively, 
51% of the shower energy is contained within a cylinder of radius 4 cm for the CORSIKA 
showers compared to 35% from the SAUND parameterisation. These fractions are approxi- 
mately independent of energy. Hence, in acoustic detectors a harder frequency spectrum for the 
acoustic signals is predicted by CORSIKA than by the SAUND parameterisation. Note that in 
the fit described in Section 5.1 the values of the parameter P 1R (equivalent to R M in equation 
20) were strongly depth dependent and much lower than the Moliere radius in water, assumed 
by the SAUND collaboration. In addition, the value of P 2 r (equivalent to s in equation 20) 
while relatively constant tended to be at a higher value (~ 1.9) than that assumed by SAUND. 



5.3 The parameterisation used by Niess and Bertin 

Hadron showers, generated by Geant4 (version 4.06 p03), were studied up to energies of 10 5 
GeV and electromagnetic showers to higher energies by Niess and Bertin [15,16]. The hadronic 
showers were parameterised as follows. 

d 2 E 

rf(z)g(r,z) (21) 



drdz 
with 

= * ( teV --e*p- te ' (22) 
A 1 \a) 

where E is the energy of the hadron shower, X is the radiation length in water, z' = z/X Q , 
b = 0.56 as determined from the fit and a is chosen to satisfy z' max = (a — l)/b. Here z' max is 
the depth in radiation lengths at which the shower maximum occurs. This is parameterised as 

z' max = 0.65 logA + 3.93 (23) 

with E c = 0.05427 GeV. The radial distribution function is parameterised as 

g(r,z)=g (^j (24) 
9 



where r,; = 3.5 cm, n = ri\ = 1.66 — 0.29(z/z max ) for r < r t and n = n 2 = 2.7 for r > r,;. 
The constant g is chosen to be (2 — rii) (n 2 — 2)/((n 2 — ni)rf ) so that the integral of the radial 
distribution is unity. 

Figure 10 shows the radial distributions from this parameterisation compared with the pre- 
dictions of CORSIKA. There is quite good agreement between the two. There is a difference in 
the normalisation with depth since Geant4, on which this parameterisation is based, produces 
showers which tend to develop more slowly with depth than those from CORSIKA (see Fig- 
ure 2). Furthermore, both this and the SAUND parameterisation (Section 5.2) assume a linear 
variation of the shower peak depth with log E whereas CORSIKA gives a clear parabolic shape 
(see equation 6). This is illustrated in Figure 1 1 . The Niess-Bertin parameterisation predicts that 
56% of the shower energy is contained within a cylinder of radius 4 cm in reasonable agreement 
with the value of 51% from CORSIKA (these values are almost independent of energy). 



6 The acoustic signals from the showers. 

The pressure, P, from a hadron shower depositing total energy E at time t resulting from the 
deposition of relative energy density e = (1/ E)(l/2-Kr)d 2 E jdrdz at a point distant d from the 
volume, dV, follows the form [13] 

p ^=£kJ e ^( s{t - d/c) ) dV <25) 



where the integral is over the total volume of the shower. Here 8 = 2.0 ■ 10 is the thermal 



expansion coefficient of the medium at 14°C, C p = 3.8 • 10 3 J kg 1 K 1 is the specific heat 



capacity and c = 1500 ms 1 is the velocity of sound in the sea water. 

Acoustic signals seen by an observer at distance r from the shower centre are computed from 
equation (25) as follows. Points are produced randomly throughout the volume of the shower 
with density proportional to the deposited energy density and the time of flight from every 
produced point to the observer calculated. The flight times to the observer are histogrammed 
over 2" bins (in this case n = 10 is chosen) centred on the mean flight time and with a suitable 
bin width, r (chosen here to be l//s). The counts in each bin of the histogram are divided by r 
yielding the function E xyz (i). The Fourier transform of the pressure wave is then 

1 f°° E 8 d 1 E 8 f 00 1 E 8 

(26) 

using the standard Fourier transform theorem, that taking the derivative in the time domain is the 
same as multiplying by iu in the frequency domain. The Fourier transform E xyz (u) at angular 
frequency co is evaluated numerically by a fast Fourier Transform (FFT) from the histogram 
E xyz (t). A correction is applied for attenuation in the water by a factor A(u) = e~ a<yU ^ r where 
a(u) is the frequency dependent attenuation coefficient. The pressure as a function of time is 
then evaluated numerically by an inverse FFT using frequency steps from zero to the sampling 
frequency (the inverse of the bin width r i.e. 1 MHz in this case). This gives 

1 n=511 

P ^ = Tm PMA{u n )e^ (27) 

n=-512 
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where Q = 27r/1024 radians and u n /27r = nfl/27r MHz is the nth frequency. The attenuation 
coefficient a(co) is computed either according to the formulae in [42] or using the complex 
attenuation given in [15, 16]. This method of calculation was computationally much faster than 
the evaluation of the space integral given in equation 18 of reference [13] and gave identical 
results. 

Acoustic pulses, computed with the complex attenuation described in [15, 16], using the 
parameterisations of the shower profile given above are shown in Figure 12. It can be seen that 
the parameterisation developed here gives similar results to that described in [15, 16] despite 
the fact that the latter was an extrapolation from low energy simulations. The parameterisation 
used by SAUND [39,40] gives smaller signals concentrated at somewhat lower frequencies. 

Further properties of the acoustic signals are shown in Figures 13 to 16. The pulses tend 
to be somewhat asymmetric with the asymmetry defined by |P TOaa; | — \P m in\/\Pmax \ + |-Pmm|- 
The complex nature of the attenuation enhances this asymmetry. This is most evident in the far 
field conditions e.g. at 5km where non complex attenuation would yield a totally symmetric 
pulse. Figure 13 shows the angular dependence of the peak pressure. Here the angle is that 
subtended by the acoustic detector relative to the plane, termed the median plane, through the 
shower maximum at right angles to the axis of the shower. The parameterisation derived here 
gives a somewhat narrower angular spread than the others. This could be due to the slightly 
longer showers predicted by CORSIKA than the others. Figure 13 also shows the asymmetry of 
the pulse as a function of this angle. The pulse initially becomes more symmetric moving out of 
the median plane and then the asymmetry becomes negative at larger angles. Figure 14 shows 
the decrease of the pulsed peak pressure with distance from the shower in the median plane and 
the asymmetry with distance in this plane. Figures 15 and 16 show the frequency composition 
of the pulses at different angles to the median plane at 1 km from the shower and at different 
distance in the median plane, respectively. 

7 Conclusions 

The simulation program for high energy cosmic ray air showers, CORSIKA, has been modified 
to work in a water or ice medium. This allows both hadron and neutrino showers to be generated 
in the medium over a wide range of energy (10 5 to 10 12 GeV). The properties of hadronic 
showers in water simulated by CORSIKA agree with those from other simulations to within 
10 — 20%. A similar uncertainty has been noted previously from the variations in CORSIKA 
showers in air generated by different models of the hadron interactions. However, none of 
the other available simulations for water cover the range of energies accessible to CORSIKA. 
The hadronic showers produced by neutrino interactions are shown to have similar profiles to 
proton showers which deposit the same amount of energy to that from the neutrino and which 
start at the interaction point of the neutrino. The properties of the neutrino interactions are 
described. A parameterisation of the shower profiles generated by CORSIKA is given. There is 
reasonable agreement with the parameterisation based on the Geant4 simulations at low energy 
(< 10 5 GeV) developed by Niess and Bertin. However, the agreement with the parameterisation 
used by the SAUND Collaboration, which is based on the NKG formalism, is less good. The 
position of the shower maximum, determined from the CORSIKA program, is found to vary 
quadratically with log E rather than linearly as assumed in the latter two parameterisations. 
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The acoustic signals generated by neutrino interactions using CORSIKA and by the two 
other parameterisations are described and their properties are studied. The acoustic signal is 
found to be very sensitive to the energy deposited close to the shower axis. 
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Interaction length of photons in water (CORSIKA) 



CD 
03 
5 



E 

o 



E10 



Ui 

c 
a 
_i 

c 
o 

o 

(0 

£ 2 

£10 2 



• Corsika - LPM Effect On 
o Corsika - LPM Effect Off 



Pair production length 

with LPM Effect 
(dash dotted curve) 




Total Interaction length 
including photonuclear 
interaction (solid curve) 



4- 



9/7 X0 Water (No LPM effect) 



8 9 10 11 12 
Log 10 Photon Energy/GeV 



Figure 1 : The interaction length for high energy gamma rays versus the photon energy measured 
in CORSIKA (data points with statistical errors). The dash dotted curve shows the pair produc- 
tion length computed from the LPM effect using the formulae of Migdal [19]. The solid curve 
shows the computed total interaction length, including both pair production and photonuclear 
interactions with the cross section from CORSIKA. The dashed line labelled 9/7X shows the 
expected pair production length without the LPM effect. Here X is the radiation length of the 
material. 
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Average dE/dz of 100 proton showers 



E 

<D 

CD 16 

N 

5 14 

LU 

■a 12 



10 



1 1 1 


1 1 1 


1 1 1 


1 1 1 


1 1 1 


1 1 1 


i i i 


i i i 


i i i 














Geant4 seawater 

Corsika seawater 



















_ 














if 

If 



























































_ 1 




















If. 

— * I ■ 


















= * / 
















^S: , i , , , 


, , , 


, , , i , , , 




T 1 r-i 


■ ■ ■ 


i i l 





200 400 600 800 1000 1200 1400 1600 1800 2000 

z (cm) 



Average dE/dz of 100 proton showers 




200 400 



600 800 1000 1200 1400 1600 1800 2000 

z (cm) 



Figure 2: Averaged longitudinal energy deposited per unit path length of 100 proton showers 
at energy 10 4 GeV (upper plot) and 10 5 GeV (lower plot) generated in Geant4 and CORSIKA 
versus depth in the water. 



16 







Depth = 250cm 



Depth = 450cm 



Depth = 650cm 



Depth = 850cm 



Depth = 1050cm 



Geant4 seawater 
Corsika seawater 



10' 
1 10 2 
10 



■= 10 3 



1 10" 



10 





,3 _ 



,2 _ 



1^10 3 



>10 2 


DC 10 
W 

"D 4 



1 10 3 



1 10" 



10 



10' 



^ 10 2 



10 



1 



,3 - 



Depth = 250cm 



Depth = 450cm 



Depth = 650cm 



Depth = 850cm 



Depth = 1050cm 



Geant4 seawater 
Corsika seawater 




5 10 15 20 25 30 35 40 5 10 15 20 25 30 35 40 



Figure 3: Averaged radial energy deposited per 20 g cm -2 vertical slice per unit radial distance 
for 100 proton showers at energy 10 4 GeV (left hand plots) and 10 5 GeV (right hand plots) 
generated in Geant4 and CORSIKA versus distance from the axis in the water for different 
depths of the shower. 
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Average y for different structure function extrapolations 
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Figure 4: The mean value of y as a function of energy for interactions computed according 
to the standard model with the PDFs of MRS99 [30], extrapolating x and Q 2 out of the fit 
range from x = 1CT 4 linearly on a log-log scale. The upper dashed curve shows the result of 
multiplying the PDFs by i.32 lo §( 10 " 4 /^) f or PDFs with x < 10" 4 and the lower dashed curve by 
dividing by this factor. The deviations of the dashed curve from the solid one is an indication 
of the precision of the standard model. 
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Figure 5: The fraction of events per unit y interval for different u M energies computed by inte- 
grating the expressions for the CC and NC cross sections. 
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Figure 6: The longitudinal distribution of the deposited energy for neutrino showers (solid) 
generated by the Herwig-CORSIKA package and proton showers (dashed) scaled to the same 
values of shower energy E w . The scaling factors applied to the average of the protons showers 
with energy 10 10 GeV were 1.0 and 3.0 for E w = 10 10 GeV and E w = 3 • 10 10 GeV, respec- 
tively. Those applied to proton showers with energy 10 11 GeV were 0.575, 1.0 and 1.65 for 
E w = 5.75 • 10 10 E w = 10 11 and E w = 1.65 • 10 11 GeV, respectively. 
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Figure 7: The solid curves show the averaged radial energy deposited per 20 g cm~ 2 vertical 
slice per unit radial distance for 70 neutrino showers with energy transfer E w = 3 • 10 10 GeV 
(left hand plots) and E w = 1.65 • 10 11 GeV (right hand plots). The incident neutrino energy 
was 2 • 10 11 GeV. For comparison the dashed curves show the distributions from proton showers 
scaled to these energies. In the left (right) hand plots protons of energy 10 10 GeV (10 11 GeV) 
were scaled by a factor of 3 (1.65). 
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Figure 8: The acoustic signal at a distance of 1 km from the shower axis in the median plane 
computed from the average of 100 CORSIKA showers each depositing a total energy of 10 9 
GeV in the water. The dotted, dashed and solid curves shows the signals computed from the 
deposited energies within cores of radius 1.025 g cm -2 , 2.05 g cm -2 and the whole shower 
(solid curve), respectively. It can be seen that most of the amplitude of the signal comes from 
the energy within a core of radius 2.05 g cm~ 2 . 



22 



10 s GeV protons 100 shower average, 10 11 GeV protons 100 shower average. 
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Figure 9: The radial distributions of the deposited energy at different depths from CORSIKA 
compared to the parameterisation used by the SAUND collaboration [39] for 10 6 GeV and 10 11 
GeV proton induced showers. 
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10 s GeV protons 100 shower average, 10 11 GeV protons 100 shower average. 
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Figure 10: The radial distributions of the deposited energy at different depths from CORSIKA 
compared to the parameterisation used by the Niess and Bertin [15, 16] (labelled NB parame- 
terisation) for 10 6 GeV and 10 11 GeV proton induced showers. 



24 



1000 



CM 



E 

o 



800 



700 



600 



o 

0) 
Q_ 

i_ 

(D 

O 



CO 500 



400 



300 



200 



- 

- 






s 

s 

/ 

y 

/ 


- 






s 

s 

s 








NB ^ " 






✓ 

SAUND ,'' 

/ 


m s ' 






s 


t S ' ^ — — 

— - m * 














S ■ 4 


"^CORSIKA 




































I i i i i 


I I I I I I I I I I I 


i i i i i i i i i i i i i 



10 



11 



12 



Log 10 E/GeV 



Figure 1 1 : The depth of the shower peak as a function of log 10 E from CORSIKA (black points) 
for the showers starting at z = 0. The solid curve shows the parameterisation according to 
equation 6. The dashed (dash dotted) lines show the values assumed by the SAUND (Niess- 
Bertin) Collaborations. 
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The Acoustic Pulse from a 10 11 GeV Shower at 1km Relative Power Spectrum of a Pulse from a 10 11 GeV Shower at 1km 
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Figure 12: The left hand plot shows acoustic pulses generated from the parameterisation de- 
scribed in section 5.1 labelled ACoRNE, the parameterisation from reference [15, 16] labelled 
NB and that from reference [39,40] labelled SAUND. These pulses were evaluated for a hadron 
shower from a neutrino interaction depositing hadronic energy of 10 11 GeV 1 km distant from 
an acoustic detector in a plane perpendicular to the shower axis at the shower maximum (the 
median plane). The right hand plot shows the frequency decomposition of the pulses in the left 
hand plot. 
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Figure 13: The left hand plot shows the variation of the peak pressure in the pulse with angle 
from the median plane at 1 km from the shower. The right hand plot shows the variation of 
the pulse asymmetry with this angle at the same distance. The curves were computed from the 
parameterisations labelled. 
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Figure 14: The left hand plot shows the decrease of the pulse peak pressure and the right hand 
plot the pulse asymmetry, both in the median plane, as a function of distance from the shower 
computed from the parameterisations. 
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Relative Power Spectrum of a Pulse from a 10 11 GeV Shower at 1km Relative Power Spectrum of a Pulse from a 10 11 GeV Shower at 1km 
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Figure 15: The left hand plot shows the frequency decomposition of the acoustic signal, com- 
puted from the parameterisation of the CORSIKA showers, at different angles to the median 
plane at a distance of 1km from the shower and the right hand plot shows the cumulative fre- 
quency spectrum i.e. the integral of the left hand plot. 



Relative Power Spectrum of a Pulse from a 1 GeV Shower Relative Power Spectrum of a Pulse from a 1 GeV Shower 




5 10 15 20 25 30 35 40 5 10 15 20 25 30 35 40 

Frequency (kHz) Frequency (kHz) 



Figure 16: The left hand plot shows the frequency decomposition of the acoustic signal, com- 
puted from the parameterisation of the CORSIKA showers, at different distances from the hy- 
drophone in the median plane and the right hand plot shows the cumulative frequency spectrum 
i.e. the integral of the left hand plot. 
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